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Abstract —Recently, communications systems that are both 
energy efficient and reliable are under investigation. In this 
paper, we concentrate on an energy-detection-based transmission 
scheme where a communication scenario between a transmitter 
with one antenna and a receiver with significantly many antennas 
is considered. We assume that the receiver initially calculates 
the average energy across all antennas, and then decodes the 
transmitted data by exploiting the average energy level. Then, 
we calculate the average symbol error probability by means 
of a maximum a-posteriori probability detector at the receiver. 
Following that, we provide the optimal decision regions. Fur¬ 
thermore, we develop an iterative algorithm that reaches the 
optimal constellation diagram under a given average transmit 
power constraint. Through numerical analysis, we explore the 
system performance. 

Index Terms —Massive MIMO, energy-detection-based trans¬ 
mission, power allocation, map detector 

I. Introduction 

The steadily growing wireless data traffic is a key reason be¬ 
hind the research for new transmission technologies with high 
performance gains and high energy efficiency. Among these 
technologies, massive multiple-input multiple-output (MIMO) 
systems, also called Large-Scale Antenna systems IT], made 
a dramatic strike since they are energy efficient and reliable 
||2l . The concept of Massive MIMO, initially proposed in 13, 
is based on using significantly many transmit and/or receive 
antennas, which leads to remarkable advances in spectral 
efficiency, beamforming gain and radiated energy efficiency 
m, 0. However, the large number of antennas poses a 
major challenge for massive MIMO to become a reality, 
since obtaining channel side information (CSI) gets more 
ponderous 0. Therefore, massive MIMO is applicable within 
time-division-duplex (TDD) systems, since channel reciprocity 
can be utilized 0. However, when applied in multi-cell 
systems such as cellular networks, massive MIMO within 
TDD systems leads to another problem, which is called pilot 
contamination. This is due to that the number of orthogonal 
pilot tones is limited in each cell, and these orthogonal pilot 
tones are reused across cells 0. 

To overcome the problem of pilot contamination, a consid¬ 
erable research effort has been expended, and several solutions 
have been proposed. One of these solutions is to design 
systems that do not require CSI at either the receiver or the 
transmitter. In this context, a fundamental work is presented in 
Q-cni where space-time coding over the Grassman manifold 
associated with the channel matrix is performed. However, 


these studies are based on the assumption of high signal- 
to-noise-ratio availability. In a more recent study El, the 
authors considered a noncoherent single-input multiple-output 
(SIMO) system with a large number of receive antennas, 
and proposed a simple energy-detection-based encoding and 
decoding scheme. They derived an upper bound for the average 
symbol error probability of the proposed scheme, and they 
showed through simulations that this upper bound will vanish 
exponentially with the increasing number of receive antennas. 
The authors also proposed a simple constellation design based 
on the minimum distance criterion. 

In another line of research, the general problem of max¬ 
imizing the system performance by deriving optimal design 
schemes is researched in numerous studies Il2]-IT4]. For 
example in na, the authors proposed a power allocation 
strategy that minimizes the bit error rate (BER) in MIMO 
spatial multiplexing systems. On the other hand, the authors 
in US developed an optimal power allocation strategy for 
orthogonal frequency-division multiplexing systems in order 
to minimize the cumulative BER. In all of these studies, it is 
assumed that both the transmitter and the receiver have perfect 
knowledge of CSI. 

In this paper, we consider an energy-detection-based com¬ 
munications system in which a transmitter with one antenna 
communicates with a receiver with many antennas. We inves¬ 
tigate a modulation and demodulation technique based on the 
calculation of energy across all receive antennas. Obtaining 
the average symbol error probability after a maximum a- 
posteriori probability (MAP) detector is employed, we provide 
the optimal decision regions. We develop an iterative algorithm 
that converges to the optimal constellation diagram under a 
given average transmit power constraint^ 

IT System Model 

As shown in Figure |l(a)| we consider a communications 
scenario in which one transmitter having a single antenna per¬ 
forms data transmission to one receiver that has N antennas. 
During data transmission, the input-output relation is given by 

y = \ix + z, 

where x is the data symbol sent by the transmitter and y is 
the N X 1-dimensional output vector at the receiver. Above, 

^ In this paper, we identify the optimal constellation diagram that minimizes 
the exact average symbol error probability rather than an upper bound to it 
as considered in HD 
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Fig. 1: System model. 


2 ; is the N X 1-dimensional additive noise vector. Each of its 
elements, Zn, is a zero-mean Gaussian random variable with 
variance cr^ for n € {1, • • • , N}. Meanwhile, h represents the 
N X 1-dimensional channel vector, each element of which, /i„, 
is also a Gaussian-distributed random variable but with mean 
/i and variance tr^. We also consider the following dynamics 
for hn'. I/ip = and cr^ = for a known real number 
K where 0 < K. Note that this is the Rician fading channel 
model QSl. We further assume that neither the transmitter 
nor the receiver knows the instantaneous realizations of the 
channel. However, they are aware of the system statistics such 
as and K. 

Since an energy-detection-based encoding and decoding 
technique is considered, the transmitter sends real positive 
symbols from a constellation of V = • • • , i/pm} as 

shown in Fig. |l(b)| rather than transmitting complex symbols. 
M is the constellation size and pm is the power level of the 
symbol, Xm, i-e., Xm = -y/Pm, for m € {1, • • • , M}. We 
note that 0 < pi < p2 < ■ ■ ■ < Pm < 00. Moreover, assuming 
that each symbol is sent with equal probability we impose 
the following average power constraint; 

1 ^ 

( 1 ) 

m—1 

where p is the average symbol power. Then, we dehne the 
signal-to-noise ratio as SNR = Finally, we denote the ratio 
of the power of the transmitted symbol to the average 
power by = ^. 

III. Average Symbol Error Probability 


We assume that in order to decode the transmitted symbol, 
the receiver utilizes the average received energy across all 
antennas, i.e., y = We can express the average received 
energy across all antennas when Xm is transmitted as 


y = 



^ Re(h* 2 ;) ,_ 

+2 — 


fLv\xj) f{y\x,) 



Fig. 2: Decision regions and boundaries for decoding. 


= P-m + + Wrn,2 + Wm,3 (2) 

— pm '^m ; 


where Wm represents the deviation of y from due to 
empirical averages. Therefore, we consider Wm as the noise 
term in our system. Now, invoki^ the central limit theorem 
and the law of large numbersj3 we show that when N 
is very large, Wm becomes a zero-mean Gaussian random 
variable with variance cr^ = j^^^i^Pm + ^ + ^Pm, i-e., 
Wm ~ A/'(0, cr^). Hence, we have y ~ A/’(//m, crp^)- It is clear 
that < (72 < ■ ■ ■ < ^M- As ssen in (l2]i, when we have a 
sufficiently large number of antennas at the receiver, we can 
characterize the average received energy with the transmitted 
symbol power p^, the channel parameters K and tr^, and 
the number of antennas N. Hence, assuming that the receiver 
applies a MAP detector, we have 


X = Xk where k = arg max f{y\xm}, (3) 


where x is the detector output. Given that Xm is transmitted, 
the conditional probability density function (pdf) of y is 


f{y\xm} = 


1 




: exp 




{y- Pmf 

2crp, 


(4) 


Noting that y is a positive real number, the receiver divides 
the positive real line into M non-overlapping decision regions: 
Di, • • • ,Dm, where Dm corresponds to the decision region 
of Xm- Following Q and (01, we can easily infer that the 
decision regions are dehned as Dm = for 3 

Am G [Pm,Pm+i] ^s depicted in Fig. |2 By default, we 
set Aq = 0 and Am = oo. Now, we can easily determine 
the symbol error probability. Initially, let us focus on the 
symbols that are at the ends of the constellation diagram. 
Assume that xi is transmitted. The receiver will be able to 
decode it correctly when y < Xi- Otherwise, the receiver 
output will be wrong. Having the conditional pdf, we can 
express the symbol error probability when xi is transmitted 

as Pe,i = f{y\xi}dy = Q -where Q{x) = 

—^ e~^ dt is the Q-function. On the other hand, when 

V'zTT ^ 21 

Xm is transmitted, the receiver will decode correctly only 


^Let {fci, • • • , Xn} be a sequence of independent and identically distributed 
random variables with mean fix and variance cr^. Defining X = EILi 
and employing central limit theorem along with the law of large numbers, we 

have limri-j-oo A/'(0,(T^) fl^ . 

























when Am -1 < V- Then, the symbol error probability is given 
by 


Pe,M = 1 



f{y\xM}dy 


= l-Q 




As for the symbols that are located in between xi and 
XM, we can easily see that when Xm is transmitted for 
mG - ,M— 1}, we have the symbol error probability 

Pe,m = 1 - / f{y\Xm}dy 

j Atti— 1 



Since each symbol is transmitted with probability we can 
express the average symbol error probability as 


M 




m—1 
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Given a power allocation vector p = [pi,--- ,pm], we can 
immediately notice that the objective function in Q, is 
separable into M—1 sub-functions, i.e.. 




M-1 


where 
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where pm < Am < Pm+i- We can show that the function 
/i(Am) in (|6]l is convex with respect to Am in its defined range 
[firn, y-m+i]- Hence, by taking the derivative of h{\m) with 
respect to Am and equating the derivative to zero, we can 
easily obtain the optimal Am that minimizes Hence, the 

optimal boundary defined with Am between two consecutive 
decision regions (i.e., the regions of pm and Pm+i) is given 
by 


Am — 


2 2 
Pm^m+l Pm+lO’m 

- rr2 


( 7 ) 


\j <7m0-m-ri {{Pm+l “ Pmf + ^Og 
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TV. Optimal Power Allocation 


In this section, we provide the optimal constellation diagram 
that minimizes the average symbol error probability Q under 
the average power constraint O for a given constellation size. 


M. Now, we have the optimization problem as 


P 




argmin Pe , 

p 


s.t. 


1 

M 


M 

Z T*™ 

m—1 


= p, 


(8) 


where p* is the vector that holds the power levels of the 
carriers of the optimal constellation diagram. The above op¬ 
timization problem does not hold a closed-form solution, and 
it is in general an NP-hard problem. Therefore, we resort to 
an iterative algorithm that converges to the optimal solution. 
Furthermore, regarding the convexity of Pg as a function of 
p, we provide the following result: 

Proposition 1: The average symbol error probability Pg 
given in (|5]l is convex in the space spanned by p. 

Proof: Omitted due to the space constraints. □ 

In the sequel, we develop a stepwise algorithm to handle 
the minimization problem in ®. Initially, let us express Pg as 
a function of p = [pi, • • • ,pm] and A = [Ai, • • • , Am-i] as 


PeipA) = j^ygiPm), ( 9 ) 

m—1 

where 

p(pm) = Q + Q (10) 

for TO G {2, • • • , M — 1}, 


9{Pi) = Q 9 (.Pm) = Q - 

Recall that pm is defined in (|2l). Now, let us assume that 
we are initially given a power allocation vector Then, 
we can easily obtain the optimal decision region boundaries 
in the first step by using (I?]), which we denote by = 
[aJ°\ • • • , A^^_]^]. Secondly, let us consider that we are given 
a vector of decision region boundaries A^°\ For any given A, 
the optimal solution for (|9l) can be obtained by the Lagrangian 
method. However, obtaining a closed-form solution for pm 
is a very difficult task, and not available without numerical 
analysis. Hence, we follow a different simplified approach, 
and treat each sub-function g{pm) separately. For given Ai, 
9 {pi) will be minimized when pi is set to 0, i.e., pi = 0. As 
for the function g{pm) defined in (fTOl i. we can see that g{pm) 
is convex with respect to pm in the defined region [Am-i, Am]- 
Then, g{pm) will be minimized when 



Pm — 


^m H” ^m— 


1 ^2 
-CT^. 


( 11 ) 


Subsequently, we can obtain pM by 


M-1 

PM = Mp - Pm- (12) 

m—1 

Then, using pi = 0, Pm in (fTTI) and pM in (fT^ . we 
form p*^^l in the second step. In the case of pM taking 
any value between pm-i and pm, i-e., Pm-i < Pm < Pm 
for TO G {2, • • • , M — 1}, we reorganize p^^l as p^^l = 
[pi, • • • ,Pm-i,PM,Pm,- ■ ■ ,PM-i]- Then, we continue the al- 
























Fig. 3: The optimum power allocation diagrams in a given 
constellation, which are obtained by employing Algorithm 1 
and brute-force when SNR = 0 dB, N = 500 and M = 6. 


gorithm until we reach a solution that satisfies the termination 
conditions. In order to formulate, we present our algorithm as 

In the following, we wrap up the above solution into an 
iterative algorithm: 

Algorithm 1: Optimal power allocation 
1: Set small e as a stopping criterion where 0 < e; 

2: Initialize p = [pi, • ‘ ‘ : Pm] and set pi = 0 such that the 
average power constraint in ([T]i is satisfied, and pi <P2 < 

■ ■ ■ < Pm', 

3: while True do 

4: Given p, compute A = [Aq, • • • , Am] using (|7]i; 

5: Given A, compute p* = [p\, - ■ ■ ,P*m\ by using (fTTT i 

and (fT2l i: 

6: if p*m-i < P*M < Pm for any m then 

7: Set p* = [p*,--- ,p*^_^,p*^,p*^,--- ,p*M_i] 

8: end if 

9: if Up - p*|p < e then 

10: break; 

11 : else 

12: Set p = p*; 

13: end if 

14: end while 

The above algorithm reaches the optimal solution after 6 
iterations when the constellation size is M = 4 and the number 
of receive antennas is = 500, and it requires about 15 
iterations to reach the solution when M = 6 with the same 
number of receive antennas. 

V. Numerical Results 

In this section, we present the numerical results. Throughout 
the paper, we consider the following settings and parameters 
unless specified otherwise. In each figure, except Fig. [2 we 
have plots as a pair in one of which we display the results when 
K = 50, and in the other the results are obtained regarding a 
channel when K = 0. We note that the channel when K = 
0 is considered to be the Rayleigh channel in which there 
is no dominant propagation along the line-of-sight, while the 
channel when AT = 50 has a strong line-of-sight propagation. 
Additionally, we set the stopping criterion e = 10“®. 



(a) K=50 



(b) K=0 


Fig. 4: The average symbol error probability, Pg, v.s. the 
number of receive antennas, N, for M = 10 and different 
SNR and K values. 


Setting AT = 50 in Fig. [2 we display the optimal power 
allocation ratios among the carriers of an optimal constellation 
diagram using Algorithm 1 when SNR = 0 dB. We further 
compare the results with those obtained using the optimal 
Brute-force search. We can easily see the good match between 
the results. Furthermore, we plot the average symbol error 
probability as a function of the number of receive antennas, 
N, in Fig. |2 when SNR = 0, 3, and 6 dB. Regardless of 
the line-of-sight character, there is dramatic decrease in the 
average probability of error with increasing TV. As expected, 
the decrease is higher when AT = 50 than it is when AT = 0. 
Similarly, we plot the average symbol error probability as a 
function of the number of constellation symbols, M, in Fig. 0 
for SNR = —6, —3, and 0 dB in order to see the effects of M 
when SNR is small. We can easily infer that the number of 
receive antennas has a great impact in obtaining small values 
of Pe even when SNR is small especially at lower constellation 
sizes. Finally, we plot the average symbol error probability, Pe, 
as a function of SNR in Fig. |2for different number of symbols 
in a given constellation, M. We see that with increasing M, Pe 
increases. We further note that due to a dominant line-of-sight 
effect we have better results in Fig. |6(a)| 

VI. Conclusion 

In this paper, we have investigated the optimal power allo¬ 
cation design for noncoherent energy-detection-based systems 
in which receivers are furnished with notably many antennas. 
Under an average transmit power constraint, we have attained 
the average symbol error probability of a MAP detector. 
























(a) K=50 (b) K=0 

Fig. 5: The average symbol error probability, Pf^, v.s. the number of symbols in a given constellation, M, for N = 500 and 
different SNR and K values. 




(a) K=50 (b) K=0 

Fig. 6: The average symbol error probability, P,., v.s. the signal-to-noise-ratio (SNR) for N = 500 and different M and K 
values. 


and we have identified the optimal decision regions for this 
setting. Showing that the average symbol error probability 
is convex in the space spanned by the symbols of a given 
constellation diagram, we have provided an iterative algorithm 
that converges to the optimal power allocation policy among 
the symbols of this constellation. Through numerical results, 
we have analyzed the effects of the channel parameters such 
as the line-of-sight character, the number of receive antennas, 
and the constellation size on the performance levels. 
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